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ABSTRACT 

The effect of angle of rotation on laminar natural convection inside the square cavity have been observed in this research. It was assumed 
that left and right walls heated isothermally, whereas other two walls act as adiabatic. This problem was solved by assuming 2-D and 
by Direct Numerical Simulation (DNS) method using ANSYS Fluent 16.0. A series of DNS simulation were carried out for different 
inclination angle (6 = 0°~90°) of the cavity at Ra = 10 3 & 10 4 . It was observed that at average Nusselt number increase up to some 
value of angle of inclination after that it decrease though this variation is not significant. 
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1. Introduction 

Study on natural convection heat transfer in square cavity 
have intense research interest due to its enormous engineering 
application like as building insulation, cooling devices for 
electronic instruments, nuclear reactor design and solar energy 
collectors etc. An extensive research have been conducted on this 
field by analytically, numerically and experimentally. 

Natural heat transfer characteristics inside the enclosure is 
controlled by some parameter, like as shape of the enclosure, 
aspect ratio, Rayleigh number and angle of inclination etc. 
Extensive study have been conducted to understand the effect of 
shape, aspect ratio of enclosure. Numerical investigation of the 
heat transfer characteristics of two dimensional rectangular 
enclosure with isothermal boundary condition for different 
aspect ratio was studied by Churchill [1] and De Vahl Davis [2]. 
An extensive review of natural convection inside the enclosure 
(non-square) were carried out by Das et al. [3]. The fluid flow 
due to the induce buoyancy force becomes turbulent when 
Rayleigh number exceed by some critical values [4]-[6]. Effect 
of inclination of rectangular cavity on the natural convection was 
studied by Sharif and Liu [7] and Kuyper et al. [8]. They 
conducted their simulation at high Rayleigh number. 

It is observed form the previous literature review that there 
is a lack of understanding about the effect of angle of inclination 
of square cavity at low Rayleigh number, which motivates for 
this research. For these reason, a series of two dimensional Direct 
Numerical Simulation (DNS) were carried for different angle of 
rotation (6 = 0°~90°) of square cavity at low Rayleigh number 
(Ra = 10 3 & 10 4 ) using ANSYS Fluent 16.0. 

2. Methodology 


assumed as adiabatic wall. To observe the effect of angle 
inclination on the natural convection in a rectangular square 
cavity, the gravitational direction was change during this 
simulation, instead of rotating the square cavity. For zero 
degree rotation (6 = 0°), gravity was considered acting along 
the negative y-axis. However, for any arbitrary rotation 6 , the 
component of gravitation forces were considered to be equal to 
—g sin 6 and — g cos 6 which is acting along x- axis and y-axis, 
respectively. During this simulation it was assumed that flow 
was steady, incompressible and two dimensional. Also 
Boussinesq approximation was considered during these DNS 
simulation run. 



Fig. 1 Domain with boundary condition 

The two dimensional governing equations (Continuity, 
Navier-Stokes equation and Energy equation) for these flow in 
Cartesian coordinates are written as follows: 


The physical system consider for this simulation is a two 
dimensional square enclosure with dimension (100 mm x 
100 mm), containing Newtonian fluid, sketched at Fig. 1. The 
left wall was consider as a rigid wall with constant temperature 
equal to 31 OK and right wall was consider as rigid wall with 
constant temperature equal to 300K. Other two walls were 
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where u, v are the velocity component in the x and y 
direction, respectively. Also, kinematic viscosity and thermal 
diffusivity represent by $ and a , respectively. 


The appropriate boundary conditions are: 


Left wall: u = v = 0 & T = 310 K 


Right wall: u = v = 0 & T = 300 K 
Lower and upper wall: u = v = 0 & — = 0 

dy 

For these DNS simulations ANSYS Fluent 16.0 were used. 
ANSYS Fluent use finite volume method to discretize the 
above governing equations. Pressure based solver was used for 
these DNS simulation run. Also, SIMPLE scheme for pressure 
and velocity coupling; Green-Gauss cell based for gradient, 
PRESTO for pressure gradient and QUICK scheme for both 
momentum and energy equations. 

In this paper, simulation were carried out for two different 
Rayleigh number (Ra = 10 3 & 10 4 ), whereas this Rayleigh 
number is defined by the following equation: 

Ra = ^(r„-r c )L 3 

da v 7 

Where, /? is thermal expansion coefficient, L distance 
among hot and cold wall. For these simulation, consider air was 
the working fluid. The different value of Ra for these DNS 
simulation run was achieved by changing the value of 
gravitational acceleration, g. Also, Local Nusselt number, Nu , 
and average Nusselt number, Nu av , were calculated with the 
following equations: 


,, hL 

Nu = — 

k 

(6) 

Nu av = Nudy 

(7) 


of iteration for Ra = 10 3 for zero degree angle of rotation. It is 
observed from Figure that average Nusselt number ( Nu av ) is 
changing significantly up to around 1000, after that it becomes 
almost stable with little variation and reach to a constant value 
which equal to 1.118 at around 1650 number iteration. Once the 
value of Nu av reach to the constant value, this simulation were 
run for more number of iteration to ensure the steady nature of 
the result. It was observed from simulation that value of Nu av 
reach to the constant value when the value of the residual of 
continuity equation, x- velocity, y-velocity and energy equation 
reach to around 10' 8 , 10" 8 , 10" 8 , 10" 13 , respectively. This 
procedure was followed for all the DNS runs carried out for this 
paper. 



Where, h is the convection heat transfer coefficient and k 
is the thermal conductivity. 

A non-uniform structural mesh was created with ANSYS 
mesh with higher resolution near to the wall, see Fig. 2. 
Specifically, for these mesh generation both horizontal and 
vertical wall were divided by 100 division with center biasing 
type with bias factor 10 to create fine mesh along the wall. This 
result was validated by comparing with previous date which is 
listed in 

Table 1. Converging criteria was set for these simulation 
based on average Nusselt number ( Nu av ) and simulation were 
run until this value of Nu av to become a constant. Fig. 3 depicts 
variation of average Nusselt number ( Nu av ) respect to number 



Fig. 3 Variation of average Nusselt number (Nu av ) with 
number of iteration at Ra = 10 3 & 0=0 


Table 1 Comparison of present work with previous work for square cavity 



Ra = 10 3 

o 

II 

Ra = 10 5 

Ra = 10 6 

Present study 

1.1179 

2.2469 

4.5347 

8.8545 

De Vahl Davis and Jones (1983) 

1.118 

2.243 

4.519 

8.799 

A. Rincon-Casado et. al. 

1.118 

2.241 

4.522 

8.819 
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3. Results and Discussion 

Fig. 4 depicts the temperature contour of the square cavity 
for different angle of rotation (6) at Ra = 10 3 . In this figure, 
temperature contour are shown without angle of inclination of 
the cavity. The exact position of the contour of the square cavity 
will be found by rotating anti-clock wise each contour shown 
in the figure with the value of 6 is mention underneath of each 
figure. It is clearly observed form the figure that at 9 = 0° at 
the top of the square cavity have higher temperature 
distribution for wider region compare to cooler region whereas 
in the bottom cooler region is wider than high temperature 
region. This result is expected due to the buoyancy effect. Fluid 
near to the left wall will rise due to the heating which reduce 
fluid density. Similarly, fluid density increase near to the right 
wall due to the cooling which produce a downward flow. 
Temperature distribution pattern are almost similar pattern for 
6 = 0°~70°, whereas at 6 = 80°~90° have different patter 
since gravitation force direction are almost perpendicular to the 
heating wall at 6 = 80°~90°, whereas it was parallel at 6 = 
0 °. 

Fig. 5 shows temperature contour of the square cavity for 
different angle of rotation ( 6 ) at Ra = 10 4 . In this figure, 
temperature contour are also shown without angle of 
inclination of the square cavity. The exact position of the 
contour of the square cavity will be found similarly by rotating 
each contour shown in the figure with the value of 0 is mention. 


It is clearly observed from figure that higher temperature zone 
is also wider in to top of square cavity and lower temperature 
zone is dominating at the bottom due to buoyancy effect which 
mention in the previous paragraph. 

Fig. 6 depicts variation of local Nusselt number (Nu) along 
the heated wall (left side wall) at different angle of inclination 
(i 6 ) for ( a) Ra = 10 3 and (b ) Ra = 10 4 . It is clearly seen from 
Fig. 6 (a) for Ra = 10 3 at 6 = 0° that Nu have higher value at 
the bottom of the heated wall which followed by slightly 
increase up to certain height and then starts to decrease and 
finally becomes almost constant near to the top of the wall. This 
similar behavior is also observed for other angle of inclination 
except 0 = 90°. The value of Nu is slightly higher at the lower 
side for 6 = 20° which decrease at 6 = 50° compare to 6 = 
0°, however no significant difference is observed at the upper 
sider of the heated wall. For 6 = 80°, Nu significantly decrease 
at the bottom of the wall and increase at the top of the wall. In 
case of 6 = 90°, Nu is equal to one along the whole heated wall 
since gravitation force acts to the perpendicular direction to the 
heated wall. Similar behavior is also observe for Re = 10 4 for 
different angle of inclination except 6 = 90°, see Fig. 6 (b ). For 
6 = 90° at Re - 10 4 , Nu have lower value at the bottom of the 
wall and higher at the top of the wall whereas as opposite 
behavior is observed for other angle of inclination. This is also 
happen since gravitation force acting to the parallel to the 
heated wall at 6 = 0° which is changing respect to the angle of 
inclination and becomes perpendicular at 6 = 90°. 
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Fig. 4 Temperature contour of the square cavity for different angle of rotation (0) at Ra = 10 3 
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(/) # = 60° (/) <9 = 70° (k) 0 = 80° (k) 0 = 90° 


Fig. 5 Temperature contour of inside the square cavity for different angle of rotation (0) at Ra = 10 4 




Fig. 6 Local Nusselt number (Nu) variation along the left wall (heated wall) at different angle of inclination (0) for (a) Re = 10 3 

and (b) Re = 10 4 


Fig. 7 shows variation of average Nusselt number ( Nu av ) 
respect to angle of inclination (6) for Re = 10 3 and Re = 10 4 . It 
is clearly observed from figure that Nu av increase up to some 
angle of rotation than again decrease. For Re - 10 3 , see Fig. 7 
(a), Nu av increase up to 6 = 20° by around 1.15% after that is 


decrease nonlinearly up to 6 = 90° and reach to the unit value 
of average Nusselt number. When Re = 10 4 , see Fig. 7 (b), Nu av 
also increase up to 6 = 40° by around 10.63% after that is 
decrease nonlinearly up to 6 = 90°. 
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(a) (b) 

Fig. 7 Effect of angle of rotation (0) on the average Nusselt number (Nu av ) at (a) Ra = 10 3 and (b) Ra = 10 4 


4. Conclusion 

This paper presents the effect of angle of orientation on the 
laminar natural convection inside a square cavity by DNS. It is 
observed that angle of orientation has positive effect on 
Average Nusselt number up to certain value, like as up to 0 = 
20° for Re - 10 3 and 6 = 40° for Re = 10 4 . Later on it decrease 
with increasing angle of orientation. However, the effect of 
angle of orientation of cavity has not much significant effect on 
Average Nusselt number. 
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